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ABSTRACT 

We present the first clear evidence for an extended dark matter halo in the Draco dwarf spheroidal 
galaxy based on a sample of new radial velocities for 159 giant stars out to large projected radii. Using 
a two parameter family of halo models spanning a range of density profiles and velocity anisotropics, 
we are able to rule out (at about the 2.5tT confidence level) haloes in which mass follows light. The 
data strongly favor models in which the dark matter is significantly more extended than the visible 
dwarf. However, haloes with harmonic cores larger than the light distribution are also excluded. When 
combined with existing measurements of the proper motion of Draco, our data strongly suggest that 
Draco has not been tidally truncated within ^ 1 kpc. We also show that the rising velocity dispersion 
at large radii represents a serious problem for modified gravity (MOND). 

Subject headings: galaxies: individual: Draco - galaxies: kinematics and dynamics - Local Group - 
dark matter - celestial mechanics, stellar dynamics 



1. INTRODUCTION 

The central velocity dispersions of many Local Group 
dwarf spheroidal (dSph) galaxies are significantly larger 
than expected for self-gravitating systems (e.g. Mateo 
1998). Assuming virial equilibrium, the implied M/ L ra- 
tios reach as high as ~ 250, making the dSphs among 
the most dark matter dominated systems in the universe. 
Given the apparent absence of dark matter in globular 
clusters (e.g. Dirsch & Richtler 1995), dSphs are also 
the smallest dark matter dominated stellar systems in the 
universe. As such, they have emerged as crucial testing 
grounds for competing theories of dark matter. 

Despite their importance, dynamical models of dSphs 
to date have been very simple. Most analyses have re- 
lied on the use of single mass, isotropic King models, with 
their associated assumptions that mass-follows-light and 
that the stellar velocity distribution is isotropic (though 
see Pryor & Kormendy 1991 and Lokas 2001 for more gen- 
eral models). Hitherto, the validity of such assumptions 
has remained unchallenged because of the small size of the 
data sets. When only small numbers of radial velocities 
are available, there is a well-known degeneracy between 
mass and velocity anisotropy (e.g., Binney & Tremaine 
1987). An increase in the line of sight velocity disper- 
sion at large radii may by due to either (1) the presence 
of large amounts of mass at large radii, or (2) tangential 
anisotropy in the velocity distribution. The primary moti- 
vation of this Letter is to break this degeneracy for Draco 
by means of improved modelling and a larger data set with 
many more stars in the outer parts. 

2. DATA 

Observations were conducted from June 23 to 26 
2000 at the 4m William Herschel Telescope using the 
AF2/WYFF0S multifibre positioner and spectrograph. A 
total of 284 stars were observed, spanning the magnitude 
range V « 17.0 to V 19.8. Of these, 159 were Draco 



members (extending to 25') with spectra of sufficient qual- 
ity to be included in our dynamical analyses. The median 
velocity uncertainty for these 159 stars was 1.9 kms~^. 
Of our stars, 62 were previously observed by Olszewski, 
Pryor & Armandroff (1996): the agreement between our 
data and the previous data is consistent with a binary frac- 
tion of ^ 40%. Full details of the data are presented in 
Kleyna et al. (2001). 

Previous studies (e.g., Hargreaves et al. 1996) have re- 
ported weak evidence of net rotational motion in Draco. 
N-body simulations show that tidal disruption of dSphs 
leads to an apparent rotation about an axis perpendicular 
to the orbit of the dSph (Oh, Lin & Aarseth 1995; Klessen 
& Zhao 2001). In our data set, we find that Draco appears 
to rotate at 6 kms~^ at a radius of 30', with a rotation axis 
position angle (PA) of 62°. However, a Monte-Carlo anal- 
ysis shows that this rotation signal is not statistically sig- 
nificant (Kleyna et al. 2001). Further, the solar rest-frame 
proper motion of Draco (Scholz & Irwin 1994) implies that 
the axis for tidally induced apparent rotation should have a 
PA of 138'^_^3^. The 3cr disagreement with the observed ro- 
tation axis argues against Draco having been significantly 
influenced by tides, a conclusion supported by the lack of 
alignment between the PA of Draco's major axis (88 ± 3°, 
Odenkirchen et al. 2001, henceforth ODOl) and its orbit. 

From our data set, it is possible to measure the line of 
sight velocity dispersion profile for a strongly dark matter 
dominated dSph for the first time. We divide the data into 
radial bins and use Bayes' theorem to obtain the proba- 
bility distribution and uncertainties of the velocity disper- 
sion {v^y^^ in each radial bin (see Kleyna et al. 2001 
for details). The top panel of Figure 1 shows a plot of 
the radial variation of (v^)^^"^. The velocity dispersion is 
clearly flat or gently rising with increasing radius. This 
plot already rules out the best-fit mass-follows-light King 
model (ODOl), for which the dispersion falls to zero at 
f'tidai = 49.4', and should have fallen to ~ 0.59 times the 
central dispersion (or ~ 5 kms~^) by 22'. 
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3. MODELLING 

3.1. Jeans Equations 

The Jeans equations for a spherical stellar system lead 
directly to a modcl-indcpcndcnt mass estimator (sec Bin- 
ney & Tremaine 1987, Eq. 4-55 & 4-56) requiring only 
knowledge of the velocity anisotropy and the true radial 
velocity dispersion (w^). Assuming spherical symmetry, it 
is straightforward to obtain {v^) from the line of sight ve- 
locity dispersion (w^) using Abel integrals. We consider 
two dispersion profiles which bracket the true situation: 
(1) a flat profile with amplitude 9.3 kms~^ (2) a profile 
rising linearly from 8.5 kms~^ in the centre to 11.25 kms~^ 
at 22'. The light is assumed to follow a Plummer law with 
a core radius ro = 9.7', as this is an excellent fit to the 
star count data (Kleyna et al. 2001). The solid and dot- 
ted curves in the lower panel of Figure 1 show the mass 
profiles obtained based on these two dispersion profiles, as- 
suming velocity isotropy. The mass enclosed within three 
core radii ranges from 6.3 - 18.0 x IO^Mq, implying M/L 
ratios of 350 — 1000, where the V-band luminosity of Draco 
is 1.8 X lO^L© V (Irwin & Hatzidimitriou 1995, henceforth 
IH95). 
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Fig. 2. — Likelihood contours of the fit of our Draco data to 
the two-parameter a, v models of Wilkinson ot af (2001). The 
contours are at enclosed two-dimensional probabilities of 
0.68, 0.90, 0.95, and 0.997. Tiic most likely value is indicated 
by a plus sign. The top and right panels represent the proba- 
bility distributions of a and v, respectively; the median of each 
distribution is represented by a square, and the triangles show 
the Icr, 2(T and Sa limits. The 2>a limits are omitted from the 
projection of a as the upper limit is strongly determined by 
the cut-off at a = 1. A 40% binary fraction is assumed. 
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Fig. 1.- Top: Line of sight velocity dispersion as a function of 
projected radius R. Tick marks show the 1 and 2rj errorbars as- 
suming a binary fraction of 40%. Bottom: Three-dimensional 
mass profiles for Draco obtained using the Jeans equations. 



3.2. Ralo Models and Distribution Functions 

We now use a two parameter family of models which 
spans a range of possibilities for the halo of Draco. The 
halo potential is given by 

-a/2 



V'(r) = 



V'o(l+r7r^)-«/^ iia^O, 
log [1 + rVr^] if a = 0. 



(1) 



The parameter a determines the mass distribution of 
Draco's halo. If a = 1, the dark matter has the same 
density distribution as the light. In models with a = 0, 
the halo has an asymptotically flat rotation curve. Finally, 
when a — —2, the visible dwarf lies in the approximately 
harmonic core of a much larger halo. The distribution 
functions for dSph density profiles in potentials (1) are 
derived in Wilkinson et al. (2001). They contain a sec- 
ond parameter v which measures the velocity anisotropy. 
Isotropic models have i/ = and models with radial (tan- 
gential) anisotropy have u > {u < 0). We generate line 
profiles from our models by numerically integrating the 
distribution functions over all transverse velocities. The 
line profile is then convolved with a binary distribution 
and Gaussian velocity measurement uncertainty profile of 
width a to give L{v,R,(j;a,iy). The probability of ob- 
serving an ensemble of points Ri with velocity Vi and ve- 
locity uncertainty Uj is then given by the product of the 
L{vi, Ri, af, a, v) values. 

Figure 2 shows the likelihood contours for the model pa- 
rameters a, z/ assuming a 40% binary fraction. The most 
likely values arc a = —0.34 and v = —0.18, indicating a 
model with significant dark matter at large radii and a 
tangentially anisotropic velocity distribution. We rule out 
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both a mass-follows-light (a = 1) and a harmonic core 
(a = —2) model at the 2.5a confidence level. Our best 
fit model implies a mass inside 3 core radii of 8^2 x lO^M©. 
The implied V-band M/ L ratio of Draco (within 3 core 
radii, or 29') is 440 ± 240. This mass profile lies be- 
tween the two curves in Figure 1. 

All our halo models (1) have cores. Cosmological sim- 
ulations favor cusped halos, such as the Navarro-Frenk- 
White (NFW) profile. The present data are insufficient 
to discriminate strongly between cored and cusped pro- 
files for Draco's dark halo - unlike the case of the Milky 
Way galaxy (Binnoy & Evans 2001). Cosmological impli- 
cations of data sets of dSph radial velocities are discussed 
elsewhere (Wilkinson et al., in prep.). 

3.3. Perigalacticon and Tidal Cut-off 




Fig. 3. Perigalacticon of Draco's orbit versus orbital ec- 
centricity. Contours are la, 2a and 3a confidence levels. 
Overplotted solid curves show the relationship between peri- 
galacticon and eccentricity assuming the labelled values of the 
tidal radii and our best-fit model for the halo of Draco. The 
dotted curve shows the earlier King model fit to Draco. The 
solid curves in the top and right panels are the probabilities 
of e and Dp, respectively, and the dashed curve in the right 
panel shows the distribution of Dp assuming no proper motion 
information. 



3.4. Perigalacticon and Tidal Cut-off 

A dSph in orbit about the Milky Way is truncated by 
the Galactic tidal field during perigalacticon passages (Oh 
& Lin 1992). Based on SDSS images of Draco, ODOl find 
no evidence for a tidal tail of Draco beyond 50' and no 
break in its luminosity profile. Thus it seems likely that 
the actual tidal truncation radius of Draco is > 1°. 

We use the radial velocity and proper motions to gener- 
ate possible Draco orbits. The proper motions are given by 
Scholz & Irwin (1994) as: cos (5 = 0.09 ± 0.05" cental, 
Us = 0.1 ± 0.05" cent~^. The line of sight velocity (cor- 
rected for the solar motion) is —98 ± 1 kms~^, while the 



distance of Draco is 82±7 kpc (Mateo 1998). We represent 
the Galactic potential by ip{r) = Vq log[(\/o^ + r'^ + a)/?"] 
(e.g., Wilkinson & Evans 1999) where the scale length 
a ~ 170kpc and the normalization vq is obtained by tak- 
ing the circular speed to be 220 kms"-"^ at the solar radius. 
We assume log-normal errors on a with 3a range 44kpc to 
489kpc; these limits correspond to masses within lOOkpc 
of 4.6 X 1O"M0 and 1.1 x IO^^Mq, respectively. Assum- 
ing Gaussian errors on all other observed quantities, we 
generate sets of initial conditions consisting of a position 
and velocity for Draco and a scale length for the Milky 
Way. Each set is integrated in the Galactic potential for a 
few orbital periods and the perigalacticon Dp and eccen- 
tricity e of the orbit is determined. Dynamical friction is 
unimportant on such short time scales. The contours in 
Figure 3 show the 1, 2 and 3a limits on the perigalacticon 
and eccentricity of Draco's orbit. While the distribution 
of eccentricities is roughly flat in the range [0.3, 0.95], the 
distribution of perigalacticons is strongly peaked at ~ 75 
kpc. This suggests that Draco is presently at or close to 
perigalacticon. 

The tidal radius for a dSph moving in the above Galactic 
potential is given by (c.f. King 1962) 

-1 



MD(rt)Dg 
Mg(l-e) 



(1 + e) 
2ev^ 



{ipp - V'a) + 



a(a' 



+ 21)2) 



(a2 + £)2)3/2 



(2) 

where Muirt) is the mass of Draco within radius rt, Mg is 
the total mass of the Galaxy, and ipp and V'a are the values 
of the potential at perigalacticon and apogalacticon. The 
solid curves in Figure 3 show the variation of Dp with e as- 
siiming a range of values for rt . Only models with rt > 3° 
fall within the la contour of Draco's perigalacticon and 
eccentricity derived from its space velocity. 

If Draco had recently been subject to strong tidal forces, 
we would expect significant elongation of the stellar dis- 
tribution, in conflict with the roundish (e — 0.3) appear- 
ance of Draco on the sky and the absence of tidal tails 
(ODOl). Chance alignment effects (e.g. Kroupa 1997) are 
unlikely to be responsible for Draco's regular appearance. 
Since any tidal extension of Draco should be aligned with 
its orbit, elongation approximately along the line of sight 
would require Draco's orbit to be similarly aligned. How- 
ever, Monte Carlo sampling of the orbital velocities using 
the observational errorbars rules out alignment within 45° 
with 98% confidence, assuming that the prior (no proper 
motion information) probability of the line of sight angle is 
uniform. The absence of statistically significant apparent 
rotation argues further against line of sight elongation. 

4. MONO AND DRACO 

As an alternative to dark matter, Milgrom (1983) pro- 
posed a modification to Newton's law of gravity (MOND) 
at low accelerations. A number of attempts have been 
made to apply MOND to the dSphs of the Milky Way 
(Gerhard 1994, Milgrom 1995). These used only the global 
velocity dispersion to obtain a MOND M/ L ratio and were 
unable to construct an unambiguous case for or against 
MOND. Lokas (2001), using published velocities, obtained 
MOND fits to the dispersion profiles of several dSphs, 
treating the global MOND acceleration scale as a free pa- 
rameter. 
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We model the luminosity density /(r) of Draco using 
a Plummer sphere with core radius tq. The velocity 
anisotropy parameter f3 (Binney 1981) is assumed to vary 
with radius as /3(r) = {-f /2)r^ /{r"^ + r^) yielding a veloc- 
ity distribution which approaches isotropy in the centre 
and becomes increasingly radial (tangential) if 7 is positive 
(negative). With these assumptions, the Jeans equations 
can be integrated to obtain the radial velocity dispersion 
crr('') as 



quantities of dark matter. 



/(r)(r2+r2)7/2 



^ ''\rl+ry^'l{r)gM{r)dr, 

(3) 

where guir) is the MOND acceleration due to the mass 
interior to radius r. The central projected V-band lu- 
minosity density Iq is 2.2 x lO^L© kpc~^ (IH95). Since 
all the mass in the MOND picture is provided by the 
stars, the Newtonian acceleration g^{r) is given simply 
by g^ir) = —GM,,{r)/r'^, where the stellar mass M,(r) 
is obtained from /(r) and depends on the M/L ratio 
of the stars. The MOND acceleration is then obtained 
via = gn/ fJ-igM/ao), where ao = 2 x 10~^ cms~^ is 
the MOND acceleration scale, /i(x) 1 for .x 1 and 
li(x) X for X ^ 0. Our assumed ao is an upper limit 
to the likely value (Sanders & Verheijen 1998), and favors 
models with lower stellar M/L ratios. 

Our analysis assumes that Draco is an isolated MOND 
stellar system. This is justified by noting that for stars be- 
tween ~ 1.7' and ~ 20' the acceleration due to the baryonic 
mass of the Milky Way disk (~ 6 x IQ-^^M©) at the distance 
of Draco is smaller than the internal accelerations. Since 
only ~ 10% of our sample lies outside 20' the assumption 
of isolation is reasonable. Representing a slightly flattened 
light distribution by a spherical model has greater impact 
when applying MOND since all the mass is in the stars. 
However, in the case of Draco this effect is not significant 
as our spherical model underestimates the mass interior to 
radius r by < 20%, resulting in M/L ratios which are at 
most 20% overestimated. 

Prom the projected dispersions, we obtain the MOND 
line profiles (assumed Gaussian) which we convolve with 
the binary velocity distribution and the error distribution 
to obtain the observable line profiles. Using these profiles, 
we analyse the velocity data set as in Section 3.2, except 
that the model parameters are now the stellar M/L ra- 
tio and the velocity anisotropy. Figure 4 shows the con- 
tours in the space of M/L and v. Again, we find that 
some tangential anisotropy is required to reproduce the 
fiat projected velocity dispersion profile. More interest- 
ingly, the most likely M/L ratio is 19 with a 3(T lower 
limit of 8.6. The assumed value of /q implies a total lu- 
minosity of 3.2 X 1O^L0 out to 3ro, approximately the 2a 
upper limit of the V-band luminosity (IH95). Our lower 
bound is therefore a true lower limit on the MOND M/L 
ratio of Draco. 

Estimates of the M / L ratio of Draco's stellar popula- 
tion based on comparison with globular clusters have a 
3cr range of 1 — 4.3 (Parmentier & Gilmore 2001, Feltz- 
ing, Gilmore & Wyse 1999). Even our 3cr lower limit on 
the MOND M/L ratio of Draco is at variance with these 
values. Thus, even by modifying the law of gravity as pre- 
scribed by MOND, we are unable to explain the observed 
motions of Draco's stars without the presence of significant 




10 20 30 40 50 60 

M/L Probability 

Fig. 4. — Likelihood contours for Draco's M/L and the 
anisotropy parameter v assuming the validity of MOND. The 
contours are at enclosed two-dimensional probabilities of 
0.68, 0.90, 0.95, and 0.997. The most likely mass to light ratio 
is M/L = I9M0/L0, and M/L = 8.6Mq/Lq is ruled out at 
the 3cr level. 



5. DISCUSSION AND CONCLUSIONS 

This Letter has presented a new set of radial veloci- 
ties for 159 stars in the Draco dSph. Our data extend to 
~ 25' and are the first observations to probe the outer- 
most regions of a strongly dark matter dominated dSph. 
The velocity dispersion profile is flat or slowly rising at 
large radii, which provides the first clear signature of an 
extended dark matter halo in any dSph. We analyse these 
data using both the Jeans equation and halo models with 
distribution functions. Our best fit model has a total mass 
within 3ro « 29' of 8^2 x IO^Mq and a slightly tangen- 
tially anisotropic velocity distribution. We are able to rule 
out the traditional mass-foUows-light models and extended 
harmonic core models at about the 2.5cr significance level. 

From orbit integrations in a Galactic potential, we ar- 
gue that Draco is currently close to perigalacticon and its 
actual tidal radius is > 1°, several times the characteristic 
length scale of the luminous matter. Draco's lack of signif- 
icant apparent rotation and the significant misalignment 
of its major axis with its orbit argue against tidal forces 
having played a major role in its evolution. Our results 
also have consequences for the Modified Newton Dynam- 
ics (MOND): a Bayesian analysis of the observational data, 
incorporating radially varying velocity anisotropy, yields a 
3cr lower hmit of 8.6 on the MOND M/L ratio still requir- 
ing the presence of a significant amount of dark matter in 
Draco. 
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